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Abstract 

A Majorana mass term for the r neutrino would induce neutrino - antineutrino 
mixing and thereby a process which violates fermion number by two units. We 
study the possibility of distinguishing between a massive Majorana and a Dirac 
r neutrino, by measuring fermion number violating processes in a deep inelastic 
scattering experiment up tX. We show that, if the neutrino beam is obtained 
from the decay of high energetic pions, the probability of obtaining "wrong sign" 
r leptons is suppressed by a factor 0{m1^6'^ /m'jj instead of the naively expected 
suppression factor E^, where Ei, is the r neutrino energy, ruu^ and are 

the r-neutrino and muon masses, respectively, and 9 is the - z/,- mixing angle. 
If rrii,^ is of the order of 10 MeV and 8 is of the order of 0.01 — 0.04 (the present 
bounds are (rrij,^ < 35MeV, 9 < 0.04) the next round of experiments may be able 
to distinguish between Majorana and Dirac r-neutrinos. 
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Introduction. Neutrino masses and mixing angles may provide one of the simplest test 
of physics beyond the Standard Model. Since neutrinos are chargeless, their mass can 
a priori be either of Dirac or Majorana type. Experimental measurements of neutrino 
masses are, however, very hard and, in fact, the only hints of nonzero neutrino masses 
come from astrophysical processes rather than from laboratory experiments. Nonzero 
neutrino masses provide a way of resolving the discrepancy between the predicted and the 
observed neutrino flux coming from the sun. Indeed, the so called MSW mechanism 
relies not only on nonzero masses but, in analogy to what happens in the quark sector of 
the electroweak theory, on nonzero mixing between the neutrinos of the first two families. 
If this mechanism were correct, very precise predictions for the masses and Cabbibo-like 
mixing angles of the first two generations could be obtained from the combination of the 
most recent experimental data0]. 

A strong constraint on neutrino masses comes from a cosmological argument . If all 
neutrinos were stable, limits on their masses may be obtained by requiring the neutrino 
density to be below the critical density needed to close the universe, Y^i^^v^ < lOOeV. 
According to this argument heavy neutrinos with masses in the range 1 KeV - 100 MeV 
can only exist assuming a sufficiently fast decay rate. Thus, one has to assume that 
the same unknown mechanism which gives the neutrinos a mass induces suitable decay 
processes turning them unstable. Moreover, if neutrinos were Dirac particles, limits on 
their masses would be also obtained from Supernova events, rriy < 30 KeV 

The bounds on the neutrino masses derived from laboratory experiments are, instead, 
much weaker than the astrophysical bounds. For the r neutrino, for example, the present 
experimental limit is niy^ < 35 MeV[^. This bound comes from an analysis of final states 
X in the r decay, r z/^X. If the r neutrino had a mass of the order of a few MeV, 
then supernova constraints would prevent it from being a Dirac particle. If, instead, it 
is a Majorana neutrino, its mass will not be restricted by supernova constraints, but it 
will have to decay sufficiently fast in order to fulfill the cosmological bound. 

In addition, if our present understanding of nucleosynthesis in the early universe were 
correct, a more stringent limit on the r neutrino mass would be obtained. As a matter 
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of fact, independently of whether the r neutrino is a Dirac or a Majorana particles, its 
mass would be bounded to be m,^^ < 1 MeV, if its lifetime were larger than (9(100 sec). 
On the contrary, if its lifetime were shorter or of the order of 10 sec, no bound on its 
mass would appear from these considerations 0. It is easy to show that this requirement 
would not be fulfilled if the only physics beyond the standard model were the neutrino 
masses and mixing angles and hence a r neutrino mass in the MeV range would imply 
either a revision of our present understanding of nucleosynthesis, or exciting new physics 
in the TeV range. 

In the following, we shall consider r-neutrinos with Majorana masses in the MeV 
range, produced in a deep inelastic experiment. Thus, we should also assume that the 
same mechanism which gives the r-neutrino a mass makes it unstable so that cosmo- 
logical and astrophysical constraints do not apply. Furthermore, one should keep in 
mind that it is always important to crosscheck astrophysical predictions with laboratory 
experiments. 

Majorana masses and mixing. As we just mentioned, for the extent of this work 
we shall consider the neutrinos to be massive Majorana particles. Therefore, neutrino- 
antineutrino mixing becomes possible. In fact, whereas a Dirac mass term 

= rriDi^LUji + h.c. (1) 

flips the spin of the neutrinos, a Majorana mass term 

'Cm = rriMi^lCuL + h.c. (2) 

mixes neutrinos with antineutrinos, where C is the charge conjugation matrix and u 
are the left handed and right handed components of the neutrino field, respectively. This 
is possible due to the fact that a Majorana field is invariant under CP conjugation so 
that transitions between particle and antiparticle are allowed, with probability ~ itt-'m- 
If one has more than two Majorana neutrinos, e.g. and z/,-, then there is in general a 
mass matrix M, 

^m = (-j,-J)cm(;;^) (3) 
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and mass eigenstates ui, 1J2 should be introduced 

z/^ \ _ / cos6' — exp(— sin \ ( 
Vr j \ exp(i(5) sin6' cos6' ) \ ^"^ 

Note that, in the above, contrary to the 2x2 Dirac type mixing matrix[|^, in addition 
to the mixing angle 9 a phase 5 appears. 

In the case of a free neutrino beam, one can show, based on simple dimensional ar- 
guments, that neutrino - antineutrino oscillations will be suppressed in comparison with 
the usual neutrino - neutrino oscillations by a factor {m^^/Eiy)^, where is the neutrino 
energy, and hence, for highly energetic neutrinos, these oscillations will be unobservably 
small p|. However, neutrino-antineutrino oscillations may be experimentally observed 
if the neutrinos proceed from a highly energetic pion beam. The physics behind this 
observation is extremely simple and is based on the helicity properties of the neutrinos as 
observed from the laboratory frame . Indeed, if one considers the frame in which the 
pion is at rest, the neutrinos will be emitted with an energy E* = (m^ — m^)/2m7r — 30 
MeV, where and are the pion and muon masses, respectively, and with definite 
helicity. When observed from the laboratory frame, however, a fraction of them, namely 
those which were moving in the backward direction with respect to the beam in the 
rest frame, will appear to have opposite helicity and, hence, they behave like Majorana 
antineutrinos. In Ref. 0, this fraction was estimated to be of the order 




If only the /i neutrino oscillations are studied, since due to experimental limits its mass 
can not be larger than a few tens of keV, then, as observed in Ref. P] , this fraction is 
unobservably small. If the /i neutrino has a small admixture with a massive r neutrino, 
with mass in the MeV range, this process, although suppressed by mixing angles, can 
be significantly enhanced. In this letter, it is our intention to make a detailed analysis 
of deep inelastic scattering experiments, up — ^ tX, in which the u beam is obtained 
from high energetic pions. We shall show how, in the light of future neutrino oscillation 
experiments, such analysis may provide very useful information about the nature of 
massive neutrinos and their masses. 
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The Experiment. The present experimental search for i/^ - z/^ mixing IIQ] is based 
on a very simple physical principle: If a high energy beam of charged pions is produced 
in the laboratory, the 7r+ (vr^) will rapidly decay into a /i+ (/i^) and a z/^ (i^ti)- If the 
neutrino beam subsequently collides with an isoscalar target, fi~ (yU"*") particles will be 
produced, which may be easily detected experimentally. Moreover, considering that the 
electroweak eigenstate z/^ is an admixture of two massive particles ui and 1^2, with 1/2 
being predominantly a particle, there will be a nonzero probability of producing r~ 
{t'^) after the collision. The detection of r particles in the final states would provide the 
first experimental clear signature of masses and mixing angles in the neutrino sector, 
and hence of physics beyond the standard model. 

Experiments like this have been carried out at Fermilab and new experiments are 
designed at CERN and Fermilab [|ll|]. The typical experimental setup is as follows: 
Assume the pions of energy (9(100 GeV) move along the z axis in the direction of the 
target. After the pion decay pipe (of order < 1 km), the neutrinos which are emitted 
mainly in the forward direction may oscillate along an oscillation pipe of length 0{1 km) 
before they hit the target, which has a typical extension (9(1-10 m^). For neutrino masses 
in the MeV range, the oscillation is too fast to be observable, so the exact length of the 
oscillation pipe is not relevant for our analysis. According to the experimental setup 
one may introduce cuts on the angle between the pion and neutrino beam, < i)cut, 
so that the neutrino beam hits the target, and on the r lepton energy, E^. > Ecut, so 
that T lepton identification is guaranteed. Furthermore, as it is usual in deep inelastic 
scattering experiments, a cut > Ql^^ with Qcut = 0{1 GeV), is introduced to ensure 
the validity of the quark parton model involved in the analysis. 

The neutrino oscillation and mixing experiment under consideration leads in general 
to excluded areas in the 6 — Am plane, where Am = 1712 — mi is the difference between 
the mass eigenvalues of the two neutrinos. For our case (m2 — 0{1 MeV), mi ~ 0) the 
oscillations are too fast, and consequently a limit 6 < 0.04 can be deduced from the non 
observation of r leptons in the Fermilab data|12|. The proposed two new experiments are 
estimated to increase the experimental accuracy by at least a factor of ten. Therefore, it 
is conceivably that r events will be observed, a small fraction of which will violate lepton 
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number conservation, so that the nature of the corresponding r neutrino mass term is 
revealed. The diagrams which describe the reactions are shown in Fig.l. Also shown in 
Fig. 1 are the corresponding reactions with a. v — v mixing Majorana propagator which 
leads to r - leptons of opposite charge. 

Amplitudes and Cross Sections. The total cross section consists on three different 
processes, the pion decay, the neutrino oscillations and finally the scattering process. 
Within a good approximation, for the case of z/^ — (z/^ — Vr) oscillations the total cross 
section factorizes, that is to say, it is a product of the pion decay rate, the oscillation 
factor and the deep inelastic lepton scattering. In the case of neutrino - antineutrino 
oscillations, the process is a little more involved and the above factorization does not 
hold. Following Ref . ||T3| , we compute the total amplitude as follows. We first consider 
the vr"*" decay amplitude to be given by 

Gp 

Mif, = -j=f^Ui{ki) [mi{l - 75) - m^(l + 75)] Vf,{k^) (6) 

were Gp and /^r are the Fermi and pion decay constants, respectively, rrii is the mass 
of the neutrino mass eigenstates z/j and /cj, and kr are the four momenta of the 
neutrino z/j, the // and r leptons, respectively. The form of the vr/iz/ interaction vertex is 
essentially dictated by the pseudoscalar nature of the pion. Analogously, for the 7r~, its 
decay amplitude is given by 

Gp 

Mi^ = --j=f^Ui{ki) [mi{l + 75) - m/,(l - 75)] v^{kfj,). (7) 

The amplitude for the scattering process of the neutrino against an isoscalar target, with 
a T~ in the final state, reads, 

Gp 

Tu,-N ^ ^ '^'^ ^ u{kr)Y{^ - lb)Ui{ki) (8) 

where is the hadronic charged current amplitude. For an antineutrino, the corre- 
sponding scattering process leads to 

T,^-N = ^<Ju> v{k,)r{l - l,)v,{h) (9) 

The total process amplitude is a product of these two amplitudes multiplied by the fac- 
tors associated to the mixing angles and the neutrino time evolution factors. Therefore, 



COS 6* — sin6'exp(— ^5) 
sinO exp{i6) cos 9 
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the total amplitude for 7r+ —>■ /x+r X is given by 

T,+^,- = Yl M,;r,^NUi^Ul expi-zEit), (10) 

i 

where Uij are the elements of the z/^ - mixing mass matrix, 

U - 

For the lepton number violating case, instead, the total amplitude reads 

r,+^,+ = J2 M.^n^M^Uir exp{-zE,t). (12) 

i 

The generalization for the vr" case is straightforward. 

From the above expressions it is easy to get the total cross sections for the processes 
by computing the square of the probability amplitude. In order to do this the hadronic 
matrix elements of the charged currents should be given. These matrix elements may 
be parametrized by |]14 



X Spins 

. l,.-^^)[,,-^^)^wr>i.,Q\ (13, 

where q = k^, — k^-, p is the nucleon momenta and u = pq/M, with M the nucleon 
mass. We treat the W exchange in the Fermi limit, i.e. we assume that = —q^ is 
much smaller than m^.This assumption is already implicit in eqs. 6 to 9. We have only 
considered the dominant structure functions, which are related to the scaling functions 

Fi{x) by 

Fi(x) = 2MWi{u,Q^) F2{x) = vW^iy^'') F^{x) = uWs{u,Q^), (14) 

where x = Q"^ /2Mv. In the parton model, the structure functions Fi{x) may be ex- 
pressed as functions of the quark densities. In the following, in order to make quantita- 



tive estimates we have used the parametrization given in ref. [|15|] . 
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One important consequence of the hadronic current structure is that the neutrino 
scattering process is enhanced in comparison to the antineutrino one. In fact, this 
can be observed, for example, in the approximation in which the charged lepton and 
neutrino masses are neglected, and ignoring small scaling violation terms appearing 
in the antineutrino cross section. The differential cross section for neutrino scattering 
against an isoscalar target is then given by 

W = (?(^) + - yf) , (15) 

while for an antineutrino it reads, 

where q{x) = x\N^{x) + iVd(x)], q{x) = x\N^{x) + iV"j(x)], and A^„(s) and N^^g^ are the 
up and down quark (antiquark) densities of the proton, respectively, and y = v/Ey 
is the fraction of neutrino energy transfered to the charged lepton. It is then obvious 
to observe, that since q{x) ^ q{x), the total neutrino cross section is enhanced by 
approximately a factor three with respect to the antineutrino one. 

Apart from the hadronic charged current matrix elements the neutrino oscillation 
process is important. For the lepton number conserving processes (I.e.) the relevant 
mass and mixing angle dependent factor is given by 

Ksl = E Ui,KU;^Ujr exp [i (Ej - E,) t] (17) 
which can be approximately expressed as 



2 cos^ e sin^ 9 



cos 



(18) 



2E, ^ 

where the second term between brackets gives the usual oscillation effect. For one of the 
masses in the MeV range, the oscillations in any of the planned experiments would be 
so fast that, once the average over the pion decay position is performed, the oscillation 
term gives no contribution to the total cross section. 

For the lepton number violating process (l.v.), instead, the oscillation term depends 
quadratically on the neutrino masses. It is given by 

K-:c = J2UiU,^U:,U,rmmjexp[t{E,-E,)t] 
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sin^ cos^ 6 



11 ( Am^ 

m2 + nil ~ 2mim2 cos ( „ „ t 



V 2^. 
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(19) 



Observe that, in this case the cross section has a nontrivial dependence on the CP 
violating parameter 5. This is a very relevant property of Majorana neutrinos, for 
which, in contrast to the Dirac neutrino cases, there is no need of mixing of the three 
generations to obtain CP violating effects. Unfortunately, it can be easily seen that, 
for neutrino mass values consistent with the present experimental constraints, the CP 
violating contribution is unobservably small and only the first term m| in the above 
expression gives a relevant contribution [nii ^ ^i)- 

Results. The total process can now be straightforwardly computed. We shall not give 
the details here, but just refer to the relevant properties. The ratio of the lepton number 
conserving to the lepton number violating process rate, in the case of initial tt^ state is 
given by 



where the first factor comes from the ratio of the oscillation factors, and the matrix 
elements Vaa' and Caa' are given by 



+ 2m,,Tr 



Ka' = 2 (m^ + ml) Tr [(1 T 75) ^.7°'" 

[IT 16) hi" ^r7"'] (21) 

Caa' = 2ml - ^J) [[(1 ± ^5) hr ^r^'] , (22) 
and an integration over phase space is understood. In order to get an idea of the order 
of magnitude of the above ratio, Eq. (|20|) , we can use the approximation k,y ~ k^. In this 
case, it is easy to show that, approximately {mi — m^^) 



p.v. ^2^ [ml + 2m2) 



(23) 



2ml {^l-ml) cr'^"^ 

Observe that, as we anticipated, the lepton number violating process is suppressed with 
respect to the lepton number conserving one by a factor of order of the ratio of the 
massive neutrino mass to the muon mass squared. Indeed, if the above approximation 
were correct, and based on the approximate magnitude of the u and u scattering cross 
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section we would get approximately a 6 % effect in the case of having a tt^ in the initial 
state, and a 0.7 % effect in the initial vr"*" case, when choosing m,^^ = 10 MeV. The actual 
numbers for the above ratio, Eq.(pOD, as we shall show below, are somewhat larger. 

In Figures 2, 3 and 4 we present the resulting ratios of cross section for different 
experimental cuts and initial pion energies. We observe that the percentage of lepton 
number violating events is roughly 1 % for tt^N and 8% for vr^A^ scattering, rather 
independent of the value of E^, (see Fig.2) and -Erc„t (see Fig. 3). From Fig. 2 it also 
follows that the ratio of the lepton number violating to the lepton number conserving 
process increases slightly if one relaxes the cut on the angle between the pion and 
neutrino beams, "^cut, while Fig. 4 shows that it depends somewhat stronger on the 
value of Ql^t. 

Conclusions. In this letter we have worked out the following idea: If the neutrinos 
have masses, the r-neutrino is probably the one with the largest mass. A priori, this 
can be either of Dirac or of Majorana type. We have discussed an experiment which 
distinguishes between a Majorana and a Dirac mass term. The main problem is that, if 
obtained from r decays, r neutrinos are high energetic and a neutrino mass determina- 
tion with high energy neutrinos is difficult since mass effects are in general suppressed 
by powers of m1^/El. Therefore, we chose the process vrA^ rX, where instead of a 
suppression proportional to the inversed r-neutrino energy squared, we found a suppres- 
sion proportional to the inverse of the muon mass squared. This is much less severe and 
allows to determine Majorana neutrino masses in the MeV range. The price we have to 
pay for this advantage is an assumption about z/^ — z/^ mixing which is needed for the 
process to take place, so that a suppression factor 9^ appears throughout in the cross 
section. For neutrino masses of the order of a few MeV, the results look promising. The 
lepton number violation process rate vr^A^ fi^r^X is of the order of a few percent of 
the lepton number conserving vr^A^ fi'^r^X' one, the exact percentage depending on 
the r-neutrino mass and the charge of the initial pion state. We have also commented 
on astrophysical and cosmological restrictions. The essential point is that cosmological 
arguments exclude stable neutrinos of mass 0{MeV). This is due to the fact that they 
would contribute too much to the present mass energy density of the universe. We have 
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ignored this bound, together with the constraints on neutrino masses coming from our 
present understanding of nucleosynthesis, by assuming that the same mechanism which 
gives the Ur a Majorana mass makes the r neutrino unstable and with a lifetime such 
that these bounds do not apply. 
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Figure Captions 



Fig. 1 

The pion nucleoli scattering processes which are discussed in this letter. Lepton 
number conserving (Dirac) processes are on the left, lepton number violating (Ma- 
jorana) processes on the right. 

Fig. 2 

The ratio of the Majorana (= lepton number violating) over the Dirac {— lepton 
number conserving) cross section as a function of the pion energy for fixed values 
of the angle cut, r energy cut and the cut in Q. On the horizontal axis the ratios 
of the Dirac cross sections for the two processes n'^N and 7r~N are given. The 

cross section for the lepton number violating process is denoted by ctm, the cross 
section for the lepton number conserving process is denoted by a^- AH results have 
a numerical error of about one percent. 

Fig. 3 

The same ratio as fig. 2, now as a function of the r energy cut. 
Fig. 4 

The same ratio as fig. 2, now as a function of the cut in Q. 
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